SRC  TR  89-81 


Analysis  Of  Coherent  Random- 
Carrier  CDMA  And  Hybrid 
WDMA/CDMA  Multiplexing  For 
High-Capacity  Optical  Networks 


By 

B.  Ghaffari  and  E.  Geraniotis 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

1989 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1989  to  00-00-1989 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Analysis  of  Coherent  Random- Carrier  CDMA  and  Hybrid 

WDMA/CDMA  Multiplexing  for  High-Capability  Optical  Networks 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Maryland, Systems  Research  Center, College 

Park, MD, 20742 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

27 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ANALYSIS  OF  COHERENT  RANDOM-CARRIER  CDMA  AND  HYBRID  WDMA/CDMA 


MULTIPLEXING  FOR  HIGH-CAPACITY  OPTICAL  NETWORKS 


B.  Ghaffari  and  E.  Geraniotis 


Department  of  Electrical  Engineering 
and  Systems  Research  Center 
University  of  Maryland 
College  Park,  MD.  20742 


ABSTRACT 

In  this  paper  we  provide  an  accurate  analysis  of  the  performance  of  a  random-carrier 
(RC)  code-division  multiple-access  (CDMA)  scheme  recently  introduced  for  use  in  high- 
capacity  optical  networks.  According  to  this  scheme  coherent  optical  techniques  are  employed 
to  exploit  the  huge  bandwidth  of  single-mode  optical  fibers  and  are  coupled  with  spread- 
spectrum  direct-sequence  modulation  in  order  to  mitigate  the  interference  from  other  signals 
due  to  the  frequency  overlap  caused  by  the  instability  of  the  carrier  frequency  of  the  laser,  or 
to  the  mistakes  in  the  frequency  coordination  and  assignment. 

The  average  bit  error  probability  of  this  multiplexing  scheme  is  evaluated  by  using  the 
characteristic  function  of  the  other-user  interference  at  the  output  of  the  matched  optical  filter. 
Both  phase  noise  and  thermal  noise  (AWGN)  are  taken  into  account  in  the  computation.  Both 
synchronous  and  asynchronous  systems  are  analyzed  in  this  context.  The  analysis  is  valid  for 
any  spreading  gain  and  any  number  of  interfering  users  and  makes  very  limited  use  of  approxi¬ 
mations.  The  performance  evaluation  of  RC  CDMA  establishes  the  potential  advantage  in 
employing  hybrids  of  WDMA  (wavelength-division  multiple-access)  and  CDMA  multiplexing 
to  combat  inter-carrier  interference  in  dense  WDMA  systems. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014-89-J-1375  and  in  part  by  the 
Systems  Research  Center  at  the  University  of  Maryland,  College  Park,  through  the  National  Science  Foundation’s  Engineering 
Research  Centers  Program:  NSF  CDR  8803012. 


I.  Introduction 


This  paper  is  motivated  by  the  recent  work  of  [1],  where  a  new  random-carrier  (RC) 
code-division  multiplexing  (CDMA)  scheme  was  introduced  for  use  in  high-capacity  optical 
networks.  According  to  this  scheme,  coherent  optical  techniques  are  employed  to  exploit 
the  huge  bandwidth  (tens  of  thousands  of  GHz)  of  single-mode  optical  fibers.  The  inherent 
instabilities  of  present-day  semiconductor  lasers  are  circumvented  by  coupling  the  optical 
multiple-access  system,  which  is  assumed  to  place  randomly  the  modulated  carriers  in 
the  available  optical  band,  with  CDMA.  In  particular,  spread- spectrum  direct-sequence 
modulation  is  employed  in  order  to  mitigate  the  interference  from  other  signals  due  to  the 
frequency  overlap  caused  by  the  instability  of  the  carrier  frequency  of  the  laser. 

The  transmitter/receiver  model  and  a  preliminary  performance  evaluation  of  this  sys¬ 
tem  were  provided  in  [1].  For  the  sake  of  simplifying  the  analysis,  several  assumptions  were 
made  in  [1]  about  the  noise  environment  (for  example,  no  thermal  noise  or  phase  noise 
were  included)  and  the  key  parameters  of  the  system  (such  as  spreading  gain,  number  of 
interfering  users)  were  assumed  infinitely  large  so  that  limiting  theorems  (e.g.,  the  Central 
Limit  Theorem)  could  be  used.  Furthermore,  the  accuracy  of  the  approximations  used 
was  not  justified  for  the  range  of  the  parameters  of  interest.  In  this  paper,  we  extend  and 
validate  the  work  of  [1]  by  providing  an  exact  evaluation  of  the  performance  of  the  RC 
CDMA  scheme  without  making  all  of  the  approximations  and  limiting  assumptions  made 
there. 

The  following  assumptions  were  made  in  [1]:  (i)  the  number  of  interfering  users  was 
taken  to  be  infinite;  (ii)  the  number  of  chips  per  bit  (spreading  gain)  of  the  spread-spectrum 
modulation  was  assumed  to  be  very  large  so  that  the  Central  Limit  Theorem  can  be 
used;  (iii)  the  error  probability  was  not  directly  evaluated  and  the  outage  probability  was 
calculated  in  terms  of  the  signal-to-noise  ratio,  which  required  the  signals  at  the  output  of 
the  optical  matched  filter  to  be  Gaussian  in  order  to  be  valid;  (iv)  the  time  delays  of  the 
various  (possibly  asynchronous)  users  did  not  enter  in  the  evaluation  of  the  performance; 
and  (iv)  the  effects  of  the  phase  noise  and  thermal  noise  on  the  system  performance  were 
not  taken  into  account.  Item  (iv)  pertains  not  only  to  the  effects  of  the  phase  noise  on  the 
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single-user  performance,  but  also  to  the  effects  of  the  phase  noise  on  the  CDMA  system 
performance  caused  by  its  effect  on  the  interference  terms  due  to  the  other  users.  Although 
all  of  the  above  assumptions  are  relaxed  in  our  paper,  in  the  calculation  of  the  characteristic 
function  of  the  multiuser  interference,  which  is  essential  for  the  approach  followed  in  our 
paper,  in  one  case  the  Central  Limit  Theorem  is  used  and  in  the  other  case  some  other 
approximation  on  the  phase  signature  sequence  is  made. 

Relaxing  the  above  assumptions  for  the  system  of  [1]  is  an  important  contribution, 
necessary  for  validating  the  usefulness  of  CDMA  in  a  variety  of  optical  network  config¬ 
urations.  Specifically,  in  some  applications  the  number  of  interfering  users  is  not  large, 
because  only  a  small  number  of  users  occuping  adjacent  frequency  bands  cause  interfer¬ 
ence.  Thus  assumption  (i)  of  the  previous  paragraph,  which  is  critical  for  the  analysis  of 
[1],  need  to  be  relaxed. 

The  analysis  presented  in  our  paper  is  valid  for  an  arbitrary  number  of  interfering 
users  and  an  arbitrary  large  number  of  chips  per  bit  of  the  direct-sequence  modulation; 
moreover,  it  takes  into  account  the  time  delays  of  the  interfering  users,  as  well  as  the 
phase  noise  of  the  lasers  (modeled  as  a  Brownian  motion  process)  and  the  thermal  noise 
present  at  the  receiver  model  (modeled  as  additive  white  Gaussian  noise  [AWGN]).  The 
performance  measure  is  the  average  probability  of  a  bit  error. 

For  the  evaluation  of  the  bit  error  probability,  we  use  the  characteristic-function 
method  intoduced  in  [2]  for  radio-frequency  (RF)  CDMA  communication  systems.  The 
evaluation  of  error  probability  in  [2]  was  carried  out  for  arbitrary  deterministic  signature 
sequences;  in  [3]  the  results  of  [2]  were  extended  to  CDMA  systems  employing  random 
signature  sequences  (i.e.,  mutually  independent  for  different  users  i.i.d  sequences  that  as¬ 
sume  the  values  +1  and  —1  with  equal  probability).  Random  signature  sequences  are  a 
very  useful  model  for  (a)  signature  sequences,  for  which  there  is  not  much  information 
available,  or  (b)  for  situations  characterized  by  a  large  number  of  potential  users  and  a 
small  number  of  users  that  are  active  or  are  expected  to  cause  interference.  In  that  last 
case,  one  can  not  assign  orthogonal  sequences  to  all  users  and  thus  needs  a  large  num¬ 
ber  of  sequences  with  reasonable  crosscorrelation  properties.  Random  sequences  make  the 
performance  evaluation  of  CDMA  systems  with  very  large  populations  feasible.  In  this 
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paper,  we  analyze  optical  RC  CDMA  sytems  with  random  signature  sequences.  In  our 
analysis,  BPSK  modulation,  as  well  as  On-Off  keying  modulation  is  used  to  modulate  the 
data  bit  stream,  while  MPSK  modulation  is  employed  for  the  signature  sequence  stream. 
Electro-optical  phase  modulators  [4]  make  these  features  feasible. 

We  evaluate  the  average  bit  error  probability  by  averaging  over  the  data  streams, 
signature  sequences, carrier  frequencies,  phases  of  the  interfering  users,  and  time  delays  (for 
asynchronous  systems).  This  is  accomplished  by  computing  the  characteristic  function  of 
the  interference  due  to  the  other  users  at  the  output  of  the  optical  filter  matched  to  a 
particular  signal.  The  computational  complexity  remains  linear  in  the  number  of  users 
and  in  the  number  of  chips  per  bit  (the  spreading  gain),  as  it  was  the  case  in  [2]  and  [3]. 
The  accuracy  of  this  computational  technique  can  be  completely  controlled  by  the  user 
and  is  determined  by  the  accuracy  of  the  integration  routines  invoked;  as  the  the  number 
of  points  in  the  integration  rule  increases,  the  required  computer  CPU  time  increases.  Any 
desirable  accuracy  can  be  attained  via  this  technique. 

Complete  numerical  results  describing  the  performance  of  the  RC  CDMA  system  will 
be  provided  in  the  near  future  and  will  be  reported  in  an  expanded  version  of  this  report. 
The  tradeoffs  between  the  various  system  parameters  (processing  gain,  signal-to-AWGN 
or  signal-to-phase  noise  ratio,  number  of  interfering  users,  and  optical  bandwidth)  will 
be  illustrated  and  interpreted.  Besides  the  average  error  probability  as  a  function  of  the 
system  parameters,  the  maximum  number  of  users  that  can  be  supported  with  this  scheme 
at  a  given  error  probability  will  be  provided.  A  comparison  with  the  results  of  [1]  will  be 
carried  out  to  examine  the  effect  of  the  approximations  and  assumptions  made  in  [1], 

As  a  byproduct  of  our  analysis,  an  exact  expression  for  the  error  probability  of  a  single 
user  coherent  optical  system  (employing  a  spreader /despreader)  disturbed  by  phase  noise 
and  AWGN  is  derived.  The  performance  of  this  system  is  also  analyzed  in  detail  in  this 
report  and  numerical  results  will  be  provided  in  the  second  version  of  the  report  for  several 
values  of  the  parameters  of  interest. 

The  performance  evaluation  of  RC  CDMA  will  validate  the  potential  advantage  in  em¬ 
ploying  hybrids  of  wavelength-division  (WDMA)  and  code-division  (CDMA)  multiplexing 
in  high  capacity  optical  networks.  A  scheme  based  on  the  principle  of  combining  the  best 
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features  of  both  FDMA  and  CDMA  multiplexing  schemes  was  proposed  recently  in  [5]. 
This  scheme  uses  primarily  WDMA  for  providing  multiple- access  capability  and  employs 
CDMA  for  protection  against  laser-frequency  instabilities  and  mistakes  in  the  frequency 
coordination  and  assignment.  This  hybrid  scheme  shows  great  potential  for  application 
to  high-capacity  optical  networks.  Our  analysis  can  be  extended  to  this  hybrid  scheme. 
Finally,  our  analysis  will  enable  us  to  compare  the  multiple-access  capability  of  this  hybrid 
WDMA/CDMA  scheme  to  that  of  a  pure  WDMA  scheme  without  CDMA  protection. 

This  paper  is  organized  as  follows.  In  section  II ,  the  model  and  the  receiver  structure 
are  described.  In  section  III ,  the  single  user  performance  for  the  BPSK  and  00K 
modulation  is  obtained.  Section  IV  extends  the  analysis  to  the  multiuser  case  in  which 
K  active  users  share  a  common  optical  channel.  The  average  bit  error  probability  for 
the  intended  user  is  obtained  by  using  the  characteristic  functions  of  the  interference  and 
AWGN .  This  section  also  contains  a  subsection  for  the  computation  of  the  characteristic 
function  of  the  other-user  interference.  Based  on  two  different  sets  of  assumptions  ,  this 
function  is  obtained  for  both  BPSK  and  00 1\  modulation.  In  section  V,  the  pdf  of  a 
useful  random  variable  which  plays  an  important  role  in  the  analysis  of  both  single  user 
and  multiuser  systems,  is  estimated.  Section  VI  extends  the  analysis  to  the  asynchronous 
system,  in  which  random  delays  are  introduced  to  the  signals  of  the  interfering  users. 

II.  Model 

There  are  I(  users  which  share  a  common  optical  channel  in  a  multiaccess  fashion. 
The  transmitted  optical  signal  is  denoted  by  S(t )  which  is  a  complex  signal  in  the  form  of 

K 

S{t)  =  M*)«m(*)ei[Wrat+* mW]  (1) 

m= 1 

where  associated  parameters  are  as  follow 

•  P  is  the  transmitted  signal  power  of  each  user. 

•  bm(t )  is  the  data  stream  of  the  m-th  user  given  by 

OO 

mo=  E 

n=  —  oo 
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where  b Lm)  denotes  the  n-th  bit  of  the  m-th  user;  6(„m)e{-l,  1}  for  BPSK  modulation  and 
&£»>e{0, 1}  for  00K  modulation.  p(t)  is  a  pulse  of  unit  amplitude  in  [0,  T], 

•  am(t )  is  the  addressing  function  or  signature  sequence  stream  used  by  user  m.  That 
is 

oo 

am(t)  =  e^m(<)  =  Y2  ti4>rnnh(t  -  nTc) 

n=— oo 

where  h(t)  is  a  pulse  of  unit  amplitude  in  [0,  Tc]  and  Tc  =  is  the  chip  duration  where  N 
is  the  number  of  chips  per  bit.  (f>mn  is  a  phase  taking  values  in  [— 7T,7r]. 

•  um  is  the  carrier  on  which  the  m-th  signal  is  sent.  This  value  is  randomly  chosen  by 
the  transmitter  laser  for  (RC)  CDMA,  or  is  preassigned  for  hybrid  WDMA/CDMA. 

•  6m(t )  is  the  phase  noise  associated  with  the  mth  transmitter  laser  which  is  a  Brownian 
motion  process  with  Lorentzian  bandwidth  /?.  The  mean  of  this  process  is  zero  and  the 
variance  is  2?r fit.  At  the  Ar-th  receiver,  the  optical  signal  S(t)  is  first  de-spread  by  a*k(t) 
which  is  the  complex  conjugate  of  ak(t)  ,  and  then  homodyne-detected  for  the  transmitted 
signal  from  user  k  (see  Fig.  1).  The  output  of  the  photodetector  is 

r(t)  =  VP  bk(t)eiA6k(t)  +  bm(t)el[w'™t+(kmit)-*kit)+A9rn{t)]  +  n(t)  (2) 


where 

i  K 

£=£ 

m  m—l 

m^k 


i  A 

W m  —  ^ m  ^ k 


A 0k(t)  =  9k(t )  -  0L(t) 


A  om(t)  =  em{t)-dL(t) 


where  dL(t)  is  the  phase  noise  of  the  local  laser  and  n(t)  is  the  complex  AWGN  process 
with  double-sided  spectral  density  ^  .  The  receiver  used  is  a  correlation  receiver  which 
is  optimum  for  the  single  user  case  with  no  phase  noise  (see  Fig.  2).  The  performance  of 
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this  suboptimum  receiver  in  the  presense  of  phase  noise  and  AWGN  is  obtained  for  the 
single  user  and  multiuser  situations  in  the  following  sections. 


III.  Single  User  Analysis 

In  this  section  the  performance  of  the  system  described  in  section  II  is  evaluated  in  terms 
of  BER  for  the  single  user  case.  The  real  part  of  the  output  of  of  the  integrator  is  denoted 
by  Y  which  is 

Y  =  b(0k)XVP  +  vVP  (3) 

where 

X  =  ^  j  cos  [A 9k(t)]dt  (4) 

and  rj  is  a  zero  mean  Gaussian  random  variable  of  variance  jpp-  The  probability  of  error 
Pe,  for  BPSK  modulation  is 


P(,  =  -Pr 

6  2 


Y  >  P\fP\b[k)  =  -1 


Y  <  p^\b[k)  =  1 


(5) 


where  p\[P  is  the  threshold.  Upon  substitution  for  Y  from  (3) 


Pe  =  \Pr[r)  >p  +  X]  +  jPr[,  <  p-X] 


(6) 


The  random  variable  X  takes  values  in  [—1, 1].  Pe  in  (6)  takes  two  possible  forms  depending 
on  the  values  of  X 


Pe  = 


|  —  |Pr  [p  —  X  <  i]  <  p  +  X/X  >  0]  w.p.  p*  =  Pr[X  >  0] 
|  +  \Pr  [p  +  X  <  Tj  <  p  —  X/X  <  0]  w.p.  q*  =  Pr[X  <  0] 


Where  ”u;.p.”  means  ’’with  probability”.  Taking  average  of  (7),  we  obtain 


Pe 


~|  ^ 

=  2  +\Pr[p  +  X<V<p-X/X<0] 
-  YPr  [p  -  X  <  n  <  p  +  X/X  >  0] 


(7) 


(8) 
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We  now  find  the  two  probabilities  in  (8).  The  first  one  is 


1_ 

q* 


Pr[p  +  X  <  r)  <  p 


X,X  <  0] 


i: 


Pr[p  +  X  <  i]  <  p  —  X]  fx(x)dx 


(9) 


where  fx(')  is  the  pdf  of  the  random  variable  X.  By  using  the  functions  $(.)  and  Q(.)  (9) 
takes  the  form 


q* 


*  (/>-*) 


2  PT 


Nn 


+  Q  \  (p  +  x) 


2  PT 


Nn 


-  1 


fx(x)dx 


(10) 


where  Q(-)  and  $(•)  are  related  to  the  standard  normal  distribution  as  follow 


<?(«)  = 

poo 

/  L_e-x2/2  •  dx 

Ja  V2n 

$(«)  = 

P  L .e-*2/2.dx 

J-O o  V2tt 

Similarly,  the  second  integral  in  (8)  is 


4  /’ 

P  J  0 


1  -  *  I  (p  -  )  -  Q  ( (p  +  X),  l2PT 


Na 


N0 


fx{x)d: 


X 


Upon  substitution  of  (10)  and  (11)  in  (8)  the  final  result  is  obtained  as 


(a) 


(12) 


The  overlines  in  (12)  indicate  expectation  with  respect  to  the  random  variable  X.  Similarly, 
for  OOK  modulation,  Pe  is  obtained  as 


Pe  = 


+  5* 


(p-X) 


(13) 
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IV.  Multiuser  Analysis 


A )  BER  Performance 

In  this  section  the  performance  of  the  system  of  section  II  is  evaluated  in  terms  of  BER 
for  the  multiuser  case.  The  output  of  the  integrator  V  is 

I 

V  =  b(0k)VPx  +  VPy}Tim  +  nVP  (14) 

m 


where 


x  =  4  f  eiAOk(t)d.t , 

T  J0 

T 

im  =  b{m)  ■  —  [ 

T  Jo 


and 


n 


=  t7fL  n(t)dt- 


(15) 

(16) 

(17) 


To  calculate  the  integral  in  (16),  we  make  some  assumptions.  First  assume  the  phase 
addressing  functions  and  <f>k(t)  in  the  chip  interval  ((n  —  1)  Tc,  nTc]  take  the  values 

4>mn  and  4>kn  in  [— 7r,7r],  respectively.  Second,  for  large  values  of  N,  it  is  reasonable  to 
assume 

A 0m(t)  =  A Om{nTc)  =  9mn  (n  -  1)TC  <  t  <  nTc  (18) 

where  6mn  is  a  zero  mean  Gaussian  random  variable  of  variance  47r /3nTc  .  Under  these 
assumptions 


b(m) 

u0 

t-m  —  jy  SZTIC 


f  ^  m  \  ei[<l>mn-<t>kn+8mn+u'm(n-l/‘2)Tc] 

'  n=  1 


(19) 


By  taking  the  real  part  of  V  in  (14) 


(20) 
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where  X  and  rj  are  defined  in  (3)  and  (4)  ,  and  Im  is 


hi?  • — - 


N 


-smc 


u'm  Tc 


N 


Y  c°s(x  mn ) 


(21) 


n—  1 


where 

Xmn  =<  (j)mn  —  4>kn  +  @mn  +  (^'rn{n  ~  l/2)Tc  >  .  (22) 

In  (22)  <  •  >  represents  [•]  mod.2ir  .  To  evaluate  the  bit  error  probability,  let  rewrite  (20) 
as 

Y  =  (b(0k)  ■  X  +  I  +  77)  s/P  (23) 

where 

/ 

i = E  *»•  (24) 

m 


The  probability  of  error  Pe,  for  BPSK  modulation,  can  be  obtained  from  (5)as  follows.  If 
we  replace  “77”  with  “I  +  t?”  equations  (6), (7)  and  (8)  are  still  valid  for  multiuser  case.  The 
first  probability  in  (8)  is 


1 

q* 


Pr[p  +  X  <  I  +  77  <  p 


X,X  <  0] 


0 

9  J- 1 


Pr  [p  +  X  <  I  +  77  <  p  —  X]  fx(x)dx 


(25) 


Upon  substitution  for  Pr  [•]  in  (25)  by  an  integral, 


q* 


J  dxfx(x)  ■  J 


p—x 

p+X 


//+„(>;)  dy, 


(26) 


where  fj+nf)  is  the  pdf  of  I  +  rj.  Since  the  pdf  //+,,(•)  is  a  real  function,  it  is  easy  to  show 
that 

fi+M  =  _  /  [Re  {$/+»,(«)}  •  cos (yu)  +  Im  {(f)i+ri(u)}  sin(yu)]  du  (27) 

Tt  JO 
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where  Re{-}  and  Im{-}  denote  the  real  and  imaginary  parts,  respectively,  of  their  complex 
arguments.  Considering  that  I  and  rj  are  independent  and  that  $i(u)  turned  out  to  be 
real  (see  part  B  of  this  section), 

1  J00 

fi+n(y)  —  ~  $i(u)  $n(u)  cos (yu)du  (28) 

TT  Jo 


Substituting  (28)  in  (26)  and  performing  the  integration  fp,  * ,  the  first  probability  in  (8) 


is 


Trq * 


poo  ‘  j*0  " 

/  $/(u)  ^(u)  cos (pu)  /  fx(x )  sin(ua;)  dx  — 

Jo  1J-1  J  u 


Using  the  same  procedure,  the  second  probability  in  (8)  is 


7 vp* 


roo  r  />1 

I  $/(«)  $r;(u)  cos  (pu)  /  fx(x)  sin(uaj)  dx 

Jo  Jo 


du 

u 


Combining  (29)  and  (30)  in  (8),  we  obtain 


1  1  f°° 

Pe  =  - - /  $l(u)  $v(u)  COS  (pu) 

2  7T  J0 


sin(udf) 


du 


u 


where 


(29) 


(30) 


(31) 


sin(uX)  =  J  fx(x)  sin(iur)  dx. 

To  evaluate  Pe  for  OOK  modulation,  we  follow  the  same  method  as  before  and  obtain 


Pe 


1  i_  r 

2  2t r  J0 


§i(u)  $n(u)  ■ 


sin(pu)  —  sin(p  —  X)i 


■  du 


u 


(32) 


In  order  to  put  (31)  in  a  more  meaningful  format  and  also  to  facilitate  the  computations, 
we  rewrite  (31)  as 


1 

7T 


f*0 o 

/  *,(«) 

Jo 


sin  (uX) 

cos  {pu) - -du 


u 


1  /  w  r  \  1  xsin(wX) 

H - /  (1  —  $/(u))  (u)  cos(pw) - -d.u. 

V  Jo  U 


(33) 
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The  third  term  in  (33)  includes  the  contribution  of  the  other-user  interference  in  BER. 
The  first  two  terms  are  single  user  contributions.  Therefore,  by  using  (13)  in  (33) 


1 

H - /  (1  —  $/(«))  $Ju)cos(pu) 

*  Jo 


sin(uX) 


u 


du. 


(34) 


Similarly  for  OOK  modulation,  Pe  is 


i  r°° 

+  27T  J0  U  “*/(«))  $l(U) 


sin(pu)  —  sin(p  —  X)u 


du 


u 


(35) 


where 


§n(u)  =  exp(- 


No  2n 


APT 


u 


(36) 


B)  Characteristic  Function  Computation 

The  computation  of  Pe  in  (34)  and  (35)  requires  the  characteristic  function  of  the  multiuser 
interference,  A>i(u)  .  This  is  defined  as 


$/(«)  =  E[eiuI]  =  E 


6(m) 

-sinc(- 


')  ££l  COS(Xmn) 


(37) 


where  Xmn  is  defined  in  (22)  .  Here,  two  approaches  are  followed  for  obtaining  this 
characteristic  function. 

According  to  the  first  approach,  it  is  assumed  that  the  carrier  frequencies  of  the  users 
are  randomly  and  uniformly  distributed  in  a  bandwith  of  W.  The  choice  of  the  phase 
signature  sequences  distribution  is  arbitrary.  This  means  that  the  sequence  { <f>mn }  takes 
values  in  [— 7r,7r]  with  arbitrary  distribution. 

In  the  second  approach,  it  is  assumed  that  the  phase  signature  sequences  are  contin¬ 
uous  and  uniformly  distributed  in  [— 7r,7r].  This  approximates  the  case  in  which  {fmn} 
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takes  M  distinct  levels  in  [— 7r,  7r]  with  probability  jj  each,  where  the  number  of  levels  M 
is  large.  In  this  case,  the  choice  of  the  carrier  frequencies  distributions  is  arbitrary. 


1. Assumption  of  Uniform  Carriers 
Let  Xmn  in  (37)  be 

Xmn  =  <^  'Ymn  T  /^mn  >  (38) 


where 

'Ymn  —  4*mn  $kn  “I-  @mn  (39) 

Pmn  =  u'm(n  -  1/2 )TC  (40) 

{o;^  }  is  assumed  to  be  iid  and  uniformely  distributed  in  a  bandwidth  W  as  big  as  10 
THZ  .  Therefore,  {/ 3mn }  for  fixed  n  are  iid  and  uniformly  distributed  in  a  bandwidth  of 
(n  —  l/2)W/(RN),  where  R  is  the  data  rate.  For  typical  values  of  R  and  N  this  value  is  still 
very  large  and  it  is  reasonable  to  consider  {<  (3mn  >}  as  iid  random  variables  uniformly 
distributed  in  [— 7r,7r]  (see  Appendix  B).  According  to  Appendix  A  ,  { Xmn }  with  respect 
to  m  are  iid  and  uniformly  distributed  in  [ — 7r,  ir]  .  Therefore,  (37)  is  expressed  as 


m 


6(”»)  t 
siu--°pj  sinc(- 


'■'>  ELi cos(x" 


.») 


(41) 


This  is  independent  of  the  index  m  .  Hence 


$l(u) 


E 


t !$sinc(2^-)  Y^,n=l  cos(Xn) 


{K-l) 


(42) 


Assuming  Yn  =  cos(Xn),  it  is  easy  to  show  that  the  sequence  {Y^}  has  zero  mean  and  zero 

correlation.  Therefore,  this  sequence  is  a  " p  —  mixing ”  (See  [6]  and  [7])  with  p(n )  =  0  . 

N 

By  using  CLT  for  large  N,  the  term  -4—  ^  cos(X„)  in  (42)  can  be  replaced  by  £  ,  where 

n= 1 

£  is  a  Gaussian  random  variable  of  zero  mean  and  variance  | .  Therefore,  (42)  becomes 


$l(u)  =  { 


E 


IU  -t= 

e  >/5v 


sinc(^)  C 


(K-l) 


(43) 
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The  expression  above  is  the  characteristic  function  of  a  Gaussian  random  variable.  There¬ 
fore, 


$i(u) 


E 


(K-l) 


(44) 


where  the  expectation  is  with  respet  to  bo  and  uj.  For  BPSK,  boe{  —  1, 1},  consequently  we 
obtain 


$l(u)  = 


(K—l) 


(45) 


For  00K,  6oe{0, 1},  consequently  we  obtain 


$/(«) 


1  1 
2  +  2eXp 


(/f-i) 


(46) 


2. Assumption  of  Uniform  Phase  Signature  Sequence 

For  the  case  in  which  the  signature  sequence  phases  are  uniformly  distributed  in  a  set  of 
equally  spaced  discrete  levels,  if  the  number  of  levels  in  this  set  is  reasonably  large,  this 
discrete  uniform  distribution  is  approximated  with  a  continuous  unifom  phase  in  [ — 7r,  7r]  . 
Let  us  express  the  characteristic  function  of  the  interference  I  as 

*/(«)  =  E*  %  By  [e"]  (47) 

Where  I  is  given  by  (24)  and  (21).  Ejr  is  the  expectation  with  respect  to  the  N(K  —  1) 
dimensional  vector  in  {Xmn}  .  Ej^  is  the  expectation  with  respect  to  the  K  —  l  dimensional 

vector  in  Finally,  E-^r  is  the  expectation  with  respect  to  the  K  —  l  dimensional 

vector  a )'  =  (u^, . . .  ,u>^).  Since  {< f>mn }  in  (22)  are  iid  and  uniformly  distributed  in  [— 7r,  7r] 
for  all  m  and  n  ,  based  on  the  result  established  in  Appendix  A,  {Xmn}  are  also  iid  and 
uniform  in  [— 7r,  7r]  for  all  m  and  n  .  Therefore, 


$/(«) 


n 


E. 


N 

n Ex 

n=  1 


•  w  »  (m  ) 

lWbo  *sin  c 

e 


cos(Xmn) 


(48) 
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We  use  the  identity 


where  Jo(-)  is  the  Bessel  function  of  the  first  kind.  Using  (49),  the  inner  expectation  in 
(48)  becomes  Jo  (b^  •  77  •  sine  which  is  independent  of  n.  Also  using  the  fact 

that  Jo(ct)  =  Jo(— a),  we  obtain  for  BPSK 


Moreover,  assuming  that  the  sequence 
m  and  (48)  becomes 


{idm}  is  iid ,  (50)  will  be  independent  of  the  index 


$/(«)  = 


(51) 


where  the  overline  is  the  expectation  with  respect  to  the  generic  u  scattered  in  the  band¬ 
width  of  the  channel.  Similarly  ,  for  OOK  modulation,  $/(u)  is  obtained  as 


$i(u)  = 


(52) 


V.  PDF  of  X 

The  pdf  of  X  which  is  defined  in  (4),  depends  on  (3T\  where  2(3  is  the  Lorentzian  bandwidth 
of  the  Brownian  motion  process  A 0*(<)  .  In  order  to  obtain  this  pdf  through  Monte  Carlo 
simulation,  (4)  is  Avritten  as 

1  F  nn6 

x  =  /  cos  (A  9k(t))dt  (53) 

n=l  d(n— 1)6 
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where  F  is  the  number  of  divisions  of  T  into  smaller  portions  of  size  8,  where  8  —  y  .  For 
F  large  enough,  let 


A 8k(t)  =  A Qk{n8)  =  9n  (n  —  1)6  <  t  <  n8 


(54) 


Upon  substitution  of  (54),  (53)  becomes 


i  F 

X  =  —  J^cos(0n) 


(55) 


n— 1 


where  { 6n  }  is  a  sequence  of  zero  mean  Gaussian  random  variables  of  variance  An/3n8  and 
correlation  An (38  min(m,  n)  .  The  statistics  of  X  is  estimated  through  the  generation  of 
the  random  sequence  { 6n }  in  a  computer  by  using  (55).  The  accuracy  of  this  estimate  is 
limitted  by  the  statistical  fluctuations  in  this  simulation.  More  specifically,  the  accuracy 
of  the  estimate  probability  of  a  quantile  q  ,  obtained  with  M  independet  simulation  trials, 

is  approximately  •  For  more  details,  the  study  of  a  similar  problem  in  [8]  is  helpful. 

In  order  to  generate  { 9n } ,  the  iteration  below  will  provide  the  required  properties  for 
this  sequence. 


6>i  =  vx 

0n  =  9n- 1  +un  2  <  n  <  F 


(56) 


where  { vn  }  is  a  sequence  of  iid  zero  mean  Gaussian  random  variables  of  varience  An [38.  It 
is  easy  to  show  that  the  sequence  { 9n }  has  the  required  properties. 


VI.  Asynchronous  System 

In  this  section  of  this  paper,  the  model  and  the  analysis  which  was  developed  for  the 
synchronous  system  in  the  previous  sections,  is  extended  to  the  asynchronous  case.  The 
transmitted  optical  signal;  S(t),  is  given  in  (1),  where  t  is  replaced  by  t—rm.  The  time  delay 
rm  is  considered  to  be  a  uniformly  distributed  random  variable  in  [0,  T]  which  represents 
the  m-th  user’s  time  delay.  At  the  A:-t,h  receiver,  the  matched  filter  is  synchronized  with  the 
A;-th  signal,  i.e.  Tk  =  0.  Therefore,  equation  (23)  is  still  valid  and  the  average  probability  of 
error  is  given  in  (34)  and  (35)  for  BPSK  and  OOK  modulation,  respectively.  Evaluation 
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of  the  characteristic  function  of  the  multiuser  interference  for  the  asynchronous  case  follows 


next.  Let  us  write  r,„  as 


where 


The  counterpart  of  (16)  becomes 


rm  =  lmTc  +  r'm 


Zm  =  L^J 


i  rT 

—  j  bm(t  Tffi  )e*[Wm  U>m  Tm  —  rm)  —  <t>k  (^)  + 

T  Jo 


The  integral  in  (58)  can  be  written  as 


rT  rTm  rT 

=  + 

Jo  Jo  J  Tm 


where 


rT m  <m  +  l  r(n-l)Te+Tln  rnTc 

J 0  n=l  J (n  —  \)Tc  n=1  J(n-1)TC+T^ 

rT  ^  .{n-1)Tc+r^  ^  nTc 

Tm  n=tm+ 2  J(n-l)Tc  n=(m  +  l  ^ (n~1)Tc+Tn 


By  using  (59),  (60)  and  (61)  in  (58),  im  becomes 


iff,  /•«“• 

rn  /  ,  ® mn  I 

1  “[  L  A«- 


(n-l)Tc+r^ 


e19™"  •  dt  +  e" 


rnTc 

j  el9mn  .  (ft; 

J(n-\)Tc+r'm 


where 


9mn  T  ^m(n  — 1) 


£5nn  —  T  4>mn  fikn  T 


emn  -  •  l[i^m+i](n)  +  bim)  ■  l[em+2, N](n) 
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emn  =  E'  •  1[1,  «„](")  +  C"  *  ' 

The  counterpart  of  (21)  is  obtained  as 


i  <  n  <  j 
other 


(66) 

(67) 


j  « 

-fm  Re{^m}  ^  ^  y  COS  Xmn  T  ^rnn  ® m  COS 


(68) 


where 


n=l 


_l  A  ;  /  u>^r' 

<*m  =  rmamc 


=  (T<  -  0*nc 


-  T-m) 


o;;  rr 

-^mn  == L^mi,‘n  l)^c  H  ^  ^mTrn  T  — 1)  ^fcn  T  9mn  > 

=<  -  1/2)TC  +  -  ^Tm  +  <f>mn  -  </>kn  +  9mn  > 

At  this  point,  similar  to  the  synchronous  case,  we  make  two  assumptions  as  follows. 


(69) 

(70) 

(71) 

(72) 


1.  Assumption  of  Uniform  Carriers 

As  before,  we  observe  that,  in  this  case  { X +ra}  and  { X~n }  with  respect  to  m  are  iid  and 
uniformly  distributed  in  [ — 7T,  7r] .  Therefore, 


*/(«)  =  IIs 


("m  E*=l  COS  Xmn  +  “m  E^l  emn  COS  Xmn) 


(73) 


This  is  independent  of  index  m.  Hence 


<s>t(u)  =  ^e  e^(“+E!=ie»cosX»++a  e:=1^c°sx„)  | 


(A'-l) 


(74) 


The  sequences  {e„  cosX„  }  and  {en  cosXn  }  are  “p  —  mixing ”  sequences.  Therefore 


$/(«)  =  {E  [e^K,'++"  *»“)]  } 


(K- 1) 


(75) 
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where  t]+  and  rj~  are  zero  mean  Gaussian  random  variables  of  variance  \  for  BPSK  and  i 
for  OOK  modulation.  Moreover,  since  r/+  and  r/~  are  uncorrelated,  they  are  independent. 
Therefore 

$/(u)  =  { E  |(A_1)  (76) 

where  r)  is  a  conditionally  zero  mean  Gaussian  random  variable  of  variance  -  )2  +  (a-  )2]  | 

for  BPSK  and  j  [(«+)2  +  (a;-)2]  for  OOK  modulation.  Finally,  for  BPSK 


$/(«)  = 


(1  -  ~)2 sinc2(u(Tc  -  r)/2)  + 


(K- 1) 


(77) 

where  the  expectation  is  with  respect  to  u  and  r.  As  before  ui  is  uniformly  distributed  in 
the  bandwidth  W  and  r  is  uniformly  distributed  in  [0,  Tc].  For  OOK  modulation,  only 
4 N  should  be  replaced  by  8 N  in  (77). 


2.  Assumption  of  Uniform  Phase  Signature  Sequence 

Since  the  two  sequences  {<j>m{n- 1)}  in  (71)  and  {</>m„}  in  (72)  are  iid  and  uniformly 
distributed  in  [ — 7T,  zt]  for  all  m  and  n,  then,  the  two  sequences  {AT+„}  and  {AT“n}  are  also 
iid  and  uniformly  distributed  in  [— 7r,7r]  for  all  m  and  n.  This  validates  equation  (73)  for 
this  case.  Moreover,  the  two  sequences  {X+n}  and  {X“n}  are  independent  of  each  other. 
To  prove  this  we  need  to  show  that  the  two  random  variables  X^in  and  X~-  are  independent 
for  all  m,n,z  and  j .  For  all  cases,  except  the  case  in  which  i  =  m  and  j  =  n  —  1,  the 
iid  and  uniform  phase  sequences  {^m(n_i)}  in  (71)  and  {4>ij}  in  (72)  provide  the  proof 
in  a  straightforward  manner.  For  the  case  m  —  i  and  j  =  n  —  1,  X+n  and  X~  are  still 
independent,  because  </>&„  and  <pk(n- 1)  are  independent  and  uniform  for  this  case.  Using 
these  facts 


$/(«) = n  e,„ 

m 


f  N 

)  n  Ex+n 

K.  n—l 

eT-amemncos^mn 

■Ex- 

mn 

-i-li.  ry  COS  X 

g  T  m  mn  mn 

l 

i 

(78) 


where  the  expectation  Em  is  with  respect  to  the  all  random  variables  with  index  m.  Upon 
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substitution  from  (49) 


$/(w)  ||  Em  /  |  Jo  (^rjiame'mn'j  ^0  y  ^  C^9) 

m  f  n — 1  ) 

For  BPSK  modulation  e+n  and  e~in  belong  to  {  —  1,1}.  Therefore,  the  inner  product  is 
independent  of  n  and  (79)  becomes 


*/(») = [j»  (%<) j-  (j«™) 


(80) 


Moreover,  assumming  that  the  sequences  {a>m}  and  {rm}  are  iid,  the  expectation  above 
will  be  independent  of  the  index  m  and  finally,  the  characteristic  function  of  the  multiuser 
interference  for  BPSK  modulation  is  obtained  as 


$/(«) 


-smc 


bJT 


2  J 


sine 


u(Tc-t) 


N 


(K- 1) 


(81) 


where  the  expectation  in  (81)  is  with  respect  to  u>  and  r,  which  are  distributed  in  \—Wf  2,  W/2]| 
and  [0,TC],  respectively. 

The  evaluation  of  the  characteristic  function  of  the  multiuser  interference  for  OOK  in 
this  case,  is  more  tedious  than  for  BPSK.  We  know  that  im  is  a  uniform  discrete  random 
variable  which  takes  values  in  {0, 1, . . .  ,  N  —  1}.  By  assuming 


(79)  becomes 


vm  ± 


*/<«) = n  £”» 

m 


'  N 

-n=fm  +  2 


(82) 

(83) 

(84) 


(SB) 
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Since  cm  ,  dm  and  vm  are  independent  of  n,  then 


*/(»)  =  ! (86) 

m 


Performing  the  expectation  above  with  respect  to  £ni ,  yields 


m 


N 


N-l 

E 


cmdm 


■At— 1- 


(87) 


By  taking  expectation  with  respect  to  (b^\b and  rearranging  the  terms,  the  final 
result  is 


1  I" 


$l(u)  =  {-  +  -[R(u>,T)R(u,,Tc-r)J 


N 


+ 


4  N 


[R{u,t)  +  R{u,Tc  -  t)] 


1  -  [R(lo,  t)R(u ,  Tc  -  r)] 
[R(u,t)R(u;,Tc-t)\ 


N 


{!<- 1) 


(88) 


where 


=  yT«nc(— ). 


(89) 
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Appendix  A 


In  this  appendix  ,  we  show  that  if  a  sequence  of  iid  random  variables,  which  are  uniformly 
distributed  in  [— 7r,7r],  are  added  to  any  other  sequence  of  random  variables  with  arbitrary 
distribution,  the  resultant  sequence  mod. 2tt  is  also  iid  and  uniformly  distributed  in  [— 7r,  7t]  . 
To  show  this,  let 

An  4>n  T  An  ^  (Al) 

where  {</>„}  is  assumed  to  be  iid  and  uniform  in  [— 7r,  7t]  for  all  n,  and  {An}  is  a  sequence 
of  random  variables  with  arbitrary  distribution.  Here  we  establish  two  claims: 

Claim  1:  Xn  and  (f)n  have  the  same  distribution. 

Proof  :  It  is  a  known  fact  that 

An|given  An  =  <f>n  (A2) 


The  pdf  of  Xn  is  obtained  as 

/x„M  =  Jj*.  l*>WA.(A)dA  (A3) 

By  usin  (A2),  (A3)  becomes 

/*»(*)  =  j^f4>n(X)f\niX)dX  =  Un(X)-  (A4) 

Claim  2:  {Anj  are  independent  for  all  n. 

Proof  :  We  show  that  any  two  random  variables  Xn  and  A;  (n  ^  l )  are  independent.  To 
show  this,  it  suffices  to  show  that  for  any  function  g(-) 

£[<7(X„MX,)]  =  E[g(Xn)]E[g(  A,)]  (A5) 

Given  An  and  A /  in  (Al),  the  first  side  of  (A5)  conditioned  on  these  values  is 

E[g(<  4>n  +  An  >).g(<  <f>i  +  A;  >)| An,  A/].  (A6) 
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Since  (j>n  and  fa  are  independent,  (A6)  becomes 

E\g(<  +  An  >)|An]  •  E[g(<  4>1  +  Aj  >)|Ai] 


(A7) 


or 


£b(X„)|A„]  •  B[j(Xi)|Ai]. 


(AS) 


From  Claim  1: 


Xn  =  Xn  |  given  An, 


Xi  =  Xi|given  Ap 


(A9) 


Therefore,  (A8)  is 


E[g{Xn)}  • 


(A10) 


This  means  that  . 

Efa(X„)s(Xi) |A„,  A,]  =  £>ls(X„)]  •  E(«(Xl)].  (A  } 

Taking  expectation  respect  to  the  A.  and  A,  gives  (A5)  and  this  completes  this  proof. 

Appendix  B 

In  this  appendix,  the  distribution  of  a  random  variable  which  is  obtained  as  the  mod. 2* 
of  a  uniform  random  variable  is  derived.  Let  assume  X  to  be  a  random  variable  winch  ,s 
uniformly  distributed  in  the  bandwidth  [0,  W\.  We  are  interested  to  know  the  distribute 
of  the  random  variable  Y  which  is  defined  as 

Y  =<X  >  .  W 


The  random 


variable  7  is  expressed  in  terms  of  some  disjoint  intervals  as  follows 

X  -  2mr  if  2mr  <  X  <  2(n  +  l)x;  0  <  n  <  M 
r  =  s  x  -  2  Mir  if  2 Mr  <  X  <  W 


(B2) 


where 

M  =  [W/2tt\ 
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For  ae[0,27r],  the  event  (Y  <  a)  is  expressed  as 


M  — 1 


|J  (X  -  2n7r  <  a,  2nvr  <X  <  2(n  +  1)tt)  |J(X  -  2Mtt  <  a,  2n7r  <  X  <  W)  (B 3) 


n=0 


The  probability  of  (B3)  is  the  summation  of  the  probabilities  of  the  indivitual  events  in 
(B3).  By  using  the  conditional  probability  and  after  some  simplifications,  Fy(oc)  ;  the 
cummulative  distribution  function  of  Y,  is 


FV(a) = 


ot\W/2'K\  +  min(ot  ,  —  27t [_Wr/27rJ ) 


W 


(£4) 


For  »  1  ,  this  distribution  approaches  to  a  uniform  one. 
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